SPACING BETWEEN PHASE SHIFTS IN A SIMPLE SCATTERING PROBLEM 



STEVE ZELDITCH AND MACIEJ ZWORSKI 
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Abstract. We prove a scattering theoretical version of the Berry- Tabor conjecture: for an almost every 
surface in a class of cylindrical surfaces of revolution, the large energy limit of the pair correlation measure 
of the quantum phase shifts is Poisson, that is, it is given by the uniform measure. 



1. Introduction and statement of the result 

The Berry- Tabor conjecture || for quantum integrable systems with discrete spectra asserts that the 
spacings between normalized eigenvalues of a quantum integrable system should exhibit Poisson statistics in 
the semi-classical limit. In particular, when the eigenvalues are scaled to have unit mean level spacing, the 
distribution of their differences should be uniform. This conjecture has been verified numerically in many 
cases [Q and has been rigorously proved in an almost everywhere sense for flat 2-tori (Sarnak [Q), flat 
4-tori (Vanderkam |fj"7f), deterministically for almost all flat tori (Eskin-Margulis-Mozes ||), and for certain 
integrable quantum maps in one degree of freedom (Rudnick-Sarnak 113], Zelditch Smilansky |ll| has 
more recently posed an analogous conjecture for scattering systems with continuous spectra. Since the 
scattering matrix S(E) at energy level E is, at least heuristically, the quantization of the classical scattering 
map, he argues that when the scattering map is integrable the eigenvalues of S(E) ( known as phase shifts) 
should exhibit Poisson statistics. In particular, he proposed that the pair correlation function of scaled 
phase shifts should be uniform for surfaces of revolution with a cylindrical end (see figure 1). The purpose 
of this paper is to prove (a somewhat modified form of) this conjecture for almost all surfaces in an infinite 
dimensional family of (pairs of) such surfaces. 

To explain the modifications and state our results, we need to introduce some notation. The surfaces 
we consider are topological discs X on which S 1 acts freely except for a unique fixed point m. The metrics 
g we consider are invariant under the S action and in geodesic polar coordinates centered at m have the 
form g = dr 2 + a(r) 2 dd 2 where a(r) defines a short range cylindrical end metric (see Sect. 2). For technical 
reasons, we are only able to analyse the phase shifts at this time in the case where g has a conic singularity 
at m. 




Figure 1. A surface of revolution with a conic singularity and a cylindrical end. 



To define the pair correlation measure, we recall that at energy A 2 the scattering matrix for a surface of 
revolution with a cylindrical end is given by a diagonal (2[A] + 1) x (2[A] + 1) matrix with entries exp(27ri<5fc(A)), 
\k\ < [A] - see Sect. 3 for a detailed presentation. The phase shifts are given by Sf-(X) and are well defined 
modulo Z. The parameter k corresponds to the angular momentum or in other words to the eigenvalues of 
the Laplacian on the cross-section (the circle in our case). When |fc| is close to A we expect no scattering 
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phenomena as the classical motion is close to the bounded motion along the cross-sections (see Fig. 2). At 
the opposite extreme, when |fc|/|A| is close to 0, the classical motion is along geodesies approaching the 
singularity on the surface. 




FIGURE 2. a) Geodesies approaching the singularity; b) Geodesies close to the cross-sections. 

Since the properties of the pair correlation measure are supposed to correspond to the properties of 
smooth classical motion it is natural, at least at this early stage, to delete the angular momenta corresponding 
to the neighbourhoods of the singularities. 

Based on this discussion we define for any e > the following measure 
d f 1 1 

Pl([a,b}) = (1 _ 2f) 2A + 1 ft { (hm,k) : Z,m,fceZ, e < |Z/A| , \m/X\ <l-c, 
(2A + 1)(1 - 2e)(S l (X) - S m (X) + k) e [a, b] } . 
In other words for / G 5(M), 

f mpKdx) = 1 1 ^ J2 f (( J - 2e ^ + 2A )(^( A ) - + k » • 

J I e J + feeZ mi , ez (12) 

e<|m/A|,|IA|<l-« v ' 

Although rather cumbersome, this definition follows the standard procedure for defining pair correlation 
measures - see the references listed above. 

Our main result concerns the (modified) pair correlation function for an infinite dimensional set Q of 
2-parameter families of surfaces of revolution 

(A, g a >P) , (a, (3) G (a„ -6,a + 6)x (-5, 6) C R 2 , 

with cylindrical ends. The precise definition of Q will be given in Proposition^ in Sect. 5. The key property 
of the metrics is that the leading parts of the phase shifts, tp a ^ , of the 2-parameter families (A, g a '@) depend 
linearly on the parameters (a, /?). This feature allows us to prove: 

Theorem. Let {g a ,p} € Q. Then for almost every pair (a, (3) (in the sense of Lesbesgue measure) and for 
any sequence {A m }^ =0 satisfying 



E 



log 3 A„ 



, A r 

m=0 



< OO . 



lim P \ {f) = / f(x)dx + /(0) , / G S(R) , e > . (1.3) 



This Theorem proves the Berry- Tabor-Smilansky conjecture for phase shifts for our class of surfaces. 
The statement can only hold almost everywhere as we can produce one parameter families of surfaces for 
which the pair correlation measure is not uniform. The proof of the Theorem is based on Proposition ^ 
below which is a somewhat stronger and more precise result. 

Acknowledgments. This note originated in a discussion between U.Smilansky and the first author at the 
Newton Institute during the program on Quantum Chaos and Disordered Systems in 1997 following a talk 
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on the results of jl8) . Smilansky pointed out that the formal WKB formula for phase shifts of a surface of 
revolution with cylindrical end closely resembled the WKB formula for eigenvalues of an integrable quantum 
map on S 2 and proposed the problem of proving that the scaled phase shifts exhibit Poisson statistics. We 
would like thank U. Smilansky for bringing this to our attention. 

The first author would also like to thank the National Science Foundation for partial support under the 
grant DMS-9703775. The second author would like to thank the National Science and Engineering Research 
Council of Canada for partial support. Both authors are grateful for the hospitality of the Erwin Schrodinger 
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2. Surfaces of revolution with cylindrical ends 

We consider a class of incomplete two dimensional smooth manifolds denoted by X \ {to}, such that X 
is a topological completion of X \ {to}. We can then consider X as a manifold with a conic singularity. The 
manifold X \ {to} is globally parametrized by (0, oo) x S 1 and we put on it metrics of revolution: 

g = dr 2 + a{r) 2 d9 2 , r E (0, oo) , 6 E S 1 . (2.1) 

The metric is assumed to be a short range cylindrical end metric, that is, we require that 

\d k ; (a(r) 2 - 1) | < C k r- 2 - k , r — > oo . (2.2) 

At to we assume a conic structure: 

o(0)=0, a' (0)^0. (2.3) 

We also make a convexity assumption by demanding that 

a'(r)>0. (2.4) 

The metric can be extended to a smooth metric on X (endowed with a natural C°° structure coming from 
polar coordinates, (r, 6) ) if and only if 

o'(0) = 1, a 2p (0) = 0, p > 0, (2.5) 



see for instance |3[. We will not assume fl2.5| ) and consequently we allow bullet like surfaces shown in Fig.l. 
The classical dynamics is given by the Hamiltonian flow of the metric: 

P=\i\l=p 2 + a{ry 2 t 2 , (2.6) 

where we parametrized T*(X \ {to}) by (x, £) = (r, 6; p,t), with p and t dual to r and 9 respectively. As is 
well known this flow is completely integrable: 

{p,t} = 0, 

and t — £(dg) is called the Clairaut integral. Abstractly, do is the vector field generating the S 1 action 
on X \ {to}. As in the case of compact simple surfaces of revolution (see we have as stronger 

statement: 

Proposition 1. For (X \ {to},<?) with the metric g satisfying (2.1)-(2.4) there exist global action angle 



variables on T*(X \ {to}). 

Although it plays no part in the proof, this is worth presenting here as the global action variables are 
closely related to the asymptotics of the phase shifts. 



Proof. The moment map 



T*(X\{m})3(x,Z)^(\(;\ gx ,Z(do))E 



9x 

-If 



has the range given by the open set B = {(61,62) : | £»2 1 < 61}. For any (61,62) E B, P _1 (6i,62) consists of 
alxS 1 orbit of a single geodesic in T*(X \ {m} (the M-action corresponds to the geodesic flow and the 
S 1 -action to the 0-rotation) . 
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In the case of case of a simple surface of revolution, the global action variables, J2), are defined by 

IiO>) = [ a, a = £-dx, (2.7) 

where (71(6), 72(b)) is a global trivialization of the bundle Hi(P^ 1 (b),'Z) of the homology groups along 
the the fibers of P. When 71(6) is chosen as the orbit of the § 1 -action, then I\ = £(dg). In the case of 
non-compact surfaces discussed here the fibers are given by R x S 1 and not by S 1 x S 1 (except for the 
degenerate case of the meridians, t = where the fiber is (R \ {0}) x S 1 where corresponds to the point 



m). Consequently the integral for I2 given by (2.7) diverges (for 71 we can still take the compact orbit of 
the S 1 -action). Hence we have to normalize the integral using the fact that the surface is asymptotic to a 
cylinder with a(r) = 1. If we take 72(b) to correspond to a geodesic in P~ l (b), then outside of the turning 
point p — (or r — for the degenerate case of the meridians) it can be parametrized by r. Then £ ■ dx 
becomes pdr and we can put 

h(b) = - lim I" (b\ - 44r) 2 dr - [ R (bj - b 2 2 )Ur , (2.8) 



7r -R^oo J \ a(r) 

that is we normalize by subtracting the "free" pdr defined by p 2 + b\ = b\ . From this we find the angle 
variables as in 0,0. □ 

3. Review of scattering theory 
There are many ways of introducing the scattering matrix on a manifold of the type we consider. Since 



we only assume (2.2), X is not a 6-manifold in the sense of Melrose - see |9[], [|3|. It is a manifold with a cusp 
metric at one end and a conic metric at the other - see || . We shall not however use this point of view here. 
Instead we will proceed more classically and we will define the scattering matrix using the wave operators 
- see H for an indication of the relation between the two approaches. As in the proof of Proposition |l| we 
need a free reference problem 

X ~ R x S 1 , go = dr 2 + d9 2 . (3.1) 
On X and X$ we define the wave groups, U(t) and Uo(t): 

U{t) : C™{X)xC™{X)3{uo,u 1 )^(u(t),D t u{t)) 

where 

(D 2 - A g )u = , t*r*=o= uq , D t u\t=o= ui . 
The operators U(t) extend as a unitary group to the energy space, H(X), obtained by taking the closure of 
C^°{X) x C^°(X) with respect to the norm 

||(uo,«i)||| = ||Vuo||i» + ||«i||i». 

The definition and properties of Uo(t) are analogous. 
We then define the M0ller wave operators 

W± : H(X ) — H(X) , 

by 

W±[w]= lim U{-t)x(r)U (t)w , w&H(X ), 

t — >±oo 

where \ G C°°([0, 00); [0, 1]), x( r ) = for r < 1 and x( r ) = 1 for r > 2, and where for r > 1 we used the 
obvious idenitification of the corresponding subsets of X and Xq. In the situation we consider the existence 
of W± is quite straightforward and we choose the wave rather than the Schrodinger picture just for variety. 
The scattering operator is 

S d =W*_W+ : H{X ) — >H(X ) (3.2) 
and, as we will see below it is a unitary operator. 
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When there is no pure point spectrum then the wave operators W± are themselves unitary. In all 
situations they are partial isometries and W± — lim t ^±oo Uq(— t)x(r)U(t). The null space of W± is the span 
of the L 2 eigenfunctions of A. Under our assumptions there could only be finitely many such eigenfunctions. 

The wave operators have the intertwining properties: 



W ± [ A ° {)= A ° i)W± 



s. 





A 90 



= 0. 



Since all operators commute with the generator of the S 1 action, dg, we decompose S using the spectral 
decompositions of A So and of dg. It is easy to check that 

UiW> 

where the Schwartz kernel of dE^ is given by 

dE° x (r,e;r',8') = £ e'"^ 9 ') ( ° 2 H e *8»W(A a -» a )*Cr-r ') (A 2 _ n *)-*dX . 

Because S commutes with the generator of the free propagator, Uo(t) we obtain the scattering matrix at 
fixed energy using the above spectral decomposition: 

S = J S(\)dEl 

and then the decomposition corresponding to the eigenvalues of dg : 

From the structure of dE® it is clear that S n (X) =0 for \n\ > |A|. For |A| = |n| we follow || and put 

S n (X)= lim S n (sgn(A)r). 

T^|A| + 

We also note that S n (X) — S- n (X). 

For \n\ < \X\, S n (X) is a unitary operator, that is, it is given by multiplication by a complex number of 
unit length: 

S n (X) = e 2mS ^ (3.3) 

and the number S n (X) is the nth phase shift at energy A 2 . Another way to think about 5(A) is as a diagonal 
unitary (2n + 1) x (2n + 1) matrix, where n = [|A|]: 

-n<k,j<n 



A more "down-to-earth" definition, following the traditional way of introducing phase shifts in one dimen- 
sional scattering, is given through asymptotic expansions in (3.6) below. 

The uniform behaviour as k and A go to infinity and k <C A is a well understood semi-classical problem. 
To describe it we separate variables in the eigen-equation of the Laplacian. We remark that this procedure 
can also provide direct proofs of the general scattering theoretical statements above. 

The Laplace operator is given by 

a' (r) 1 

A g = D 2 - i—jJ-D r + —r\2 D e 
a(r) a(r) z 
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and on the eigenspaces of Dg it acts as 

a(r) ' ' a(r) 2 



„o a'(r) „ 1 o 
A„ = Dl - -i—-!-D r + — -^n 2 



, - n 2 2a"(r)a{r) - (o'(r)) . 



, , ( 3 - 4 ) 

a(r) 2 4a(r) 2 

The reduced operator appearing in brackets in the second line above has a self-adjoint realization on 
L 2 ((0, oo) r ) and for large A it can be considered semi-classically: 

A„ - A 2 = X 2 a(r)-^P(x, h)a(r)? , h = , x = ^ , 

P{x, h) = (hD r ) 2 + V(r; x, h) - 1 , (3.5) 



J 2a"(r)a(r)-( a'(r)) 2 T/ ^.^def, 
a(r) 2 4a(r) 2 



V(r;x,h) = ^—-h 2 ~" v/ " v/ . , 2 V " V// , ^o(r;x) = V(r;x,0). 



The principal symbol of P(x, h) is given by p = p 2 + x 2 /a(r) 2 — 1 and the natural range of x for which 
semi-classical methods are applicable is given by < e < \x\ < 1 — e. In fact, since a(r) is one at infinity, we 
approach zero energy when x 2 is close to 1. On the other hand when x —> the characteristic variety of p 
has a singular limit - see Fig. 3. A detailed analysis of the x — » limit has to involve the lower order terms 
in V(r;x, h). In particular, miraculous cancellations in the expansions due to the interaction between the 
leading and lower order terms occur when we have product type conic singularities since we can then use the 
theory of Bessel functions. The general situation is, at least to the authors, unclear at the moment. What 
is quite clear is that we have a uniform expansion in h/x. 




Figure 3. The characteristic variety of P(x,h). 

The phase shifts 6 n (X) are related to the semi-classical phase shifts of the operator P(x, h), ip(x, h) which 
are defined by asymptotics of solutions: 

'I s 
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We recall now the essentially standard asymptotic properties of tp - see |L0[ , Chapter 6 and for a more 
microlocal discussion JTlJ] . 

Proposition 2. As h — * 0, V'O*'; ft.) defined by (3.6) ft.as an asymptotic expansion uniform in e < \x\ < 1 — e 
/or any /lied e > : 

^(x, /i) ~ ^(s) + h*- + h 2 ip 2 (x) + ■■■ (3.7) 



2 



/■OO j 

^(^)=/ ((l-V (r, 2 ;))|-(l- a; 2 )5jdr. (3.8) 



□ 

We remark that when we translate the asymptotics to the coordinates on T*(X \ {to}): x = t/X, 
X 2 = p 2 + t 2 /a(r) 2 we obtain the second action variable defined in the proof of Proposition [j]. 

As mentioned in the introduction, we can describe this connection between the phase shifts and action 
variables by saying that S(X) is a quantum map on G(X,g), the space of geodesies. We now digress to 
explain this statement in more detail. For simplicity we will consider S at integral values of A and denote 
them by N. 

Since it is not needed in the calculation of the limit pair correlation function we give a somewhat sketchy 
discussion and refer to |Q for background on Toeplitz quantization. See also |l5| for a related discussion 
from a physicist's point of view. 

We can identify G(X,g) with the set S* n (X ro ) of incoming vectors at the parallel 

, , def ^ r f -1 

X ra = X n {r = r } . 

As in Proposition || we have to to delete e-neighborhoods of the singular set, given by {|t/A| < e} where 
t = I\ = £(dg) is the first action variable and A 2 = p 2 + i 2 a(r) _1 is the energy. 

We denote the deleted space of geodesies by G £ (X, g) and identify it with the set S* ne (X ro ) of incoming 
vectors at X ro with incoming angle satisfying \6\ > e and \9 — w\ > e. If we consider the deleted space of 
geodesies as a phase space then, on the quantum level, it corresponds to the sequence of truncated Hilbert 
spaces H.N.e spanned by the eigenfunctions {e m6 } of the quantum action with e < < 1 — e where 
l\ = —idg. Here, 1/iV plays the role of the Planck constant and we restrict it to integral values. Since Ti.N,e 
is invariant under S(N) we may restrict the latter to a unitary scattering matrix S £ (N) on H.N,e- 

We now state the somewhat informal: 

Proposition 3. The sequence {S e (N)} is a semiclassical quantum map over G e {X,g) associated to the 
classical scattering map 

P '■ ^in.e^ro) ~~ * S in e (X ro ) 

where [3(x, £) is obtained by following the geodesic J( x ,£) through (x, £) until it intersects X ro for the last time 
and reflecting the outgoing tangent vector inward. 

Proof. From the explicit formula 

S £ (N) = (e 2 ^^6 kJ U lk/Nllj/Nl < {1 ^, S k (N) = ^ _0 ( 3. 9) 

we see that S e (N) is the exponential of N times the Hamiltonian 

on TLn,h where Xe is a smooth cutoff function defining the truncated Hilbert space. 

The truncated phase space G e (X,g) is symplectically equivalent to a truncated S 2 , equipped with its 
standard area form, with neighbourhoods of the poles and of the equator deleted. Indeed, the equivalence 
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is defined by the identity map between global action-angle charts on the surfaces. This map intertwines the 
obvious S 1 actions which rotate the spaces. The quantization of this chart then defines a unitary equivalence 
on the quantum level which intertwines the operators dg on cylinder and sphere (they can be considered as 
the angular momentum operators). The equivalence is specified up to a choice of 2N + 1 phases by mapping 
the spherical harmonic of degree N which transforms under rotation by 8 on S 2 by e zke to the exponential 
e lke with k £ [— N, . . . , N]. The map is completely specified by requiring that the spherical harmonic be real 
valued along 8 = 0. 

Thus we may identify -ff/v,e with a Hamiltonian over the compact phase space S 2 . Since it is a function 
of the (Toeplitz) action operator I\ it is necessarily a semiclassical Toeplitz operator of order zero with 
principal symbol Xe(Ii I E)ip(Ii I E) on S 2 . The semiclassical parameter iV is identified in the Toeplitz theory 
with a first order positive elliptic Toeplitz operator with eigenvalue N in 7ijv ~~ see [ fl8| and references given 
there. Hence NHn,6 is a first order Toeplitz operator of real principal type. As in the essentially analogous 
case of pseudodifferential operators, the exponential of a first order Toeplitz operator of real principal type 
is a Fourier integral Toeplitz operator whose underlying classical map is the Hamilton flow generated by ip. 

We now wish to identify this map at time one with the classical scattering (or billiard ball) map on 
S* n e (X ro ), that is, we wish to prove that /3 = exp where 3^ is the Hamilton vector field of ip. 

Indeed, let us work in the symplectic action-angle coordinates (0,li) where 8 is the angle along X ro . 
The Hamilton flow of ip then takes the form 

exptS^(0,ii) = (0 + *w,Ji), u = d h ^. (3.10) 

At time t = 1 the angle along the parallel X ro changes by to. We claim that ui is also the change in angle 
along the incoming geodesic through 8 G X ro in the direction I\ as it scatters in the bullet head before 
exiting again along X ra . 

To see this, we use Proposition || which shows that ip is closely related to the second action variable: in 
the notation of the proof of Proposition [j] 

h(h,b 2 ) = b 1 r/j(^j , h(h,b 2 ) = b 2 . 

Since b\ = A is preserved by the flow we can fix it at A = 1 and then 

On the other hand the equations of motion show that 

dO 8 2ta{r)- 2 I x 



+ 



1 - 77^1 ■ ( 3 - 12 ) 



dr r 2p a(r) 2 \ a(r) 

It follows that lu(Ii) is twice the change in angle as the radial distance changes from rg to its minimum along 
the geodesic. 

The piece of the geodesic lying in the bullet-head consists of two segments: the initial segment beginning 
on S ro and ending upon its tangential intersection with the parallel X r _/j t ^ closest to m, and the segment 
beginning at this intersection and ending on X ro . The change in 0-angle along both segments is the same, 



so that the total change in angle during the scattering is given by the integral (3.12) above. This shows 
that (3 and exp 5^ have precisely the same formula in action-angle variables and completes the proof of the 
proposition. □ 
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4. Exponential sums 



Following we will reduce the study of ( Jl . 2| ) to a study of certain exponential sums. We first remark 
that because of symmetries of 5k (A) we can study a slightly simpler expression 

Pl(f) = (1 \ )x E E .f((l-e)X(S 3 (X)-6 k (X)+m)) 
as one easily checks that 



roeZ e <j/A,fc/A<l-e 



f 



P\(f) = oPM / o 



1 



The reduction to exponential sums follows from an application of the Poisson summation formula in m: 

2 

JlS k (X) 



[(l-2e)A] 2 £^((l 



2w£ 



(1 - 2e)A 



E < 

Kk/\<1- 



From this we see that 



Px(f) = /(0) + / /(2<R + OA^oc(l) + E{(f) 



em) d = f 



1 



[(l-2e)AP 



^ez\{o} 



(1 - 2e)A 



E 



e> «(« fc (A)-« J (A)) 



(4.1) 



E<fc/x,i/A<i- 



Ideally, we would like to show that E e x (f) — > as A — > 0. That however seems very hard and for some 
surfaces is simply not true. As in [[[jj , fl2|| , |is|| we will instead consider families of surfaces and resulting 
families of scattering phase shifts: 



Replacing 4(A) by S^(X) in (4.1) we define 



(ao - 7, a + 7) x (-7, 7) 3 (a, 0) ^— > 5^(X) . 
(a - 7, a Q + 7) x (-7, 7) 9 (a, f3) 1 — ► £ A (/; a, /3) 



(4.2) 



To see the point of doing this we recall from |14 ,JL2| and [p| the following simple 
Lemma 1. If for any f e 



"0+7 /-7 



I^a(/; a, [3)\ 2 dad[3 < C e jF(X) , 



1 a — 7 ^—7 

i/ien /or any sequence {A m }^ =0 smc/j f/iai 



E F ( A ^) < 

m=0 



E{ m (f;a,f3)^0, m — > 00 , V/e 
almost everywhere in (a, /3) £ (ao — 7, c*o + 7) x ( — 7, 7)- 

□ 

When Su'^(X) have a somewhat idealized form, the crucial estimate comes from |is|| , Theorem 5.1.1 
where it is loosely based on the Vinogradov method. Since we will need a further development of these 
estimates we present a slightly modified proof. 
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Proposition 4. If in (4.1) and (p 



6%' f3 (X) = ak + (3\$[ t) . *eC°°((0,l)), \&'\( e ,i-e) \>C e >0; 



then for any f G <S(M) 



Proo/. Let G C° 



|^(/;a,/?)| 2 dad/3-0/. e 



log 3 A 
A 



, A 



have the following properties 

ps(t) > supp p s C (-<5, 5) . 

The estimate of the lemma will clearly follow from 

/ / ps(a)ps((3)\El(f-a : P)\ 2 dadf3 = o( ] ^ 

Using the representation of E^, (4.1), the left hand side of ( |4.3| ) can be rewritten as 

\4 E E 



(4.3) 



A 4 



E 



E 



e<fc 1 /A,j 1 /A<l- e e<k 2 /\,j 2 /\<l- 
^1 t^JI fc 2 



ft - fci) - l 2 {j 2 - k 2 ))p s A 4 f J -f t $ T T 



A 



A 



J2 



A 



A 



(4.4) 



where we dropped the overall factor of (1 — 2e)~ 2 and put <?(£) = /(27r(l — 2e)~ 1 £). From now on we will 
drop the parameter e altogether: we can for instance extend $ as a strictly convex or concave function to 
[0, 1] adding additional positive terms to the sums which are being estimated or we can shift and rescale the 
variables. 

The support condition on ps implies that 

- h) - £ 2 (j 2 - k 2 ) = 



A 



< 



(4.5) 



A 



To understand the second expression we apply the mean value theorem twice to the difference of $'s. For 
that we make a simple observation: if <f>" has a fixed sign on [a — e, b + e] then if (m + h) /2, (m—h)/2 G [a, b] 

lmin|(/>"| dZ lmax|0"| 
2 max |0" | - dm™' ' ~ 2 mm\c/)"\ ' 
In our case we put <f> = <f>(»/A) and hi = ji — ki 7^ and m, = + fej. Then with 



we have 



and (4.5) implies 



(x,h( m ) d = ^{m, h) 



l/C < d m Z Xth (m) <C, 



hhi=l 2 h 2 , 0< \hi\ < A 



where we can invert ^A.h 2 m view of (|4.6|). 



< C<5A , < mi < 2A , 



(4.6) 



(4.7) 
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Thus we want study the sets of six integers (hi, rrii, U), i = 1,2, satisfying (4.7). We first note that, say, 
/i2 is determined by h\,l\,li. Then we see from the second inequality in (4.7) and from (4.6) that for fixed 
{h\,hi,m\,t\) there are 

A 



0(1) max 1, 



\hhi 



m2 's satisfying (4.7). 



When A > \l\h\\ the contribution to (4.4) is estimated by 



A 



^4 5Z 5Z ^2 a2 1/ /, i 

Krai<2A -*<hi<* £ 1 =^0«2#0 1 1 



A 



ff l A 



< CS 



log A 



i«7(ob / Isolde 

lil>A 141 



log 2 A 



When A < then the number of iri2 , s is uniformly bounded for each choice of the other variables. 

We want to count the triples (hi,hz,£i) satisfying the first equation of (4.7) as a function of £2 and A. Let 
F(A,^2) denote that number. If d(n) denotes the number of divisors of n 7^ then 

F(X,t 2 )<8 Y <W 2 |) 

0<h 2 <\ 

since t\h\ — £2^2 and each factorization into a product has to be counted twice since t\ and h\ can be 
interchanged. Then 

G(A,iV) d = f Y ^(V 2 )<8]T Y <W*) 
<C Y d(n) 2 <C'XN (log X + log N) 3 , 

l<n<NX 

by a theorem of Ramanujan - see 0, Sect. 18.2 and references given there. Hence the part of (4.4) corre- 
sponding to the bounded number of rn2's is bounded by 



C\- 4 Y m**\9\Y, F M') 



Km 1 <2\ 



<C'X- 3 / (G(A,A|£|)+l)|«/(« 



< C a 



log 3 A 



A 



where we used summation by parts and then approximation by the Riemann integral. This completes the 
proof of the proposition. □ 



We now recall from Proposition that for surfaces we consider we have 



1 1 



where for e < x < 1 — e, d^tp2 = Cfc, <:(!)■ 



A 



k 1 



6 k (X) = Xtp T +7 + 7^2 V'T > e< T <1 - £ ' A > 00 



A' A 



A 



(4.8) 



When the fam ily of surfaces depends on two parameters, (a, (3), so that (2.\)-(2A) hold uniformly, then 
we also have ( |4.8| ) uniformly with respect to the parameters. Hence to apply Proposition || to our case we 



need to estimate the contribution of the error terms coming from . That is given in 
Proposition 5. Let 6^ (A) be given by (4.8) uniformly in (a, (3) € (cto — 7, a>o + 7) x (—7, 7) with 



The 



i> a >P(x) = ax + (3$(x), $eC°°((0,l)), |*"r (£ ,i- e ) I >Ce >0. 



Q o+7 n /bg 3 A 



Q0-7 



|^(/;a,/3)| J dad i 9 = / , ( 
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Proof. We observe that ip a ^ is defined for all (a, 0) and that we can extend S^' 13 to all (a, 0) by smoothly 
cutting off the lower order terms for (a, 0) glxl \ (—7 — e, 7 + e) x (a — 7 — e, ao + 7 + e)- Using the 
inequality |x| 2 < 2|y| 2 + 2|x — y| 2 and Proposition || we see that we need to estimate 

2 



^0 



; «(5^(A)-5<*'"(A)) _ i£(^-"(A)-«5 j "- 3 (A)) 



i<i/l,j/A<l- 



ps(a)ps((3)dad(3 . 



where we use the notation of the proof of Proposition H. We now introduce 



and put 



t(z) = 2i sin Z - exp (-d) , e ix - e*» = e lx r(x - y) , ^ E C c 
2 V 2/ z 

A f i f , a ,0 fk 1\ ,a,0 { j 1 



€ V A 



= ^ t r 7 ,T -fi 



A A 



A' A 



From Proposition g and the obvious properties of r we see that ^,j,fe,A is C°° and that it satisfies the following 
estimates 

Pa+Pfi 



<C P ( 1 



A 



(4.9) 



Hence we have to look at 



A 2 ^— ' V A / V A 



e<fc/A,j/A<l-e 



ps(a)ps(0)dad/3 . 



We would like to proceed as in the proof of Proposition [| but now taking of Fourier transforms has to 
be replaced by integration by parts. Using the analysis of the differences of <5^ ,/3 presented there we need to 
estimate 



A 2 £ £ 



h (I 



in a 



«2 / «2 

A 5 A 



(£ 1 /A) 2 (4/A> 2 /(A^ 1 ^ 2 ), 



(4.10) 



where 



2 [A] 2 [A] [A] [A] 

I{\l U l 2 )=Y J J2 £ £ ^l^-^2^2)" 2 <A- 1 (^^2a/ l2 (™2)-^l/iia/ ll (™l))> _2 , 



mi — lm2- 1 ' 1 i = -[ A l >*2 = -[A] 



and where, as is usual, we write (x) = (1 + x 2 )^ . We note that we can absorb the terms (£i/X) 2 and {I2/X) 2 
into the g terms. 

We first observe that 



^(A-5m)- 2 = 0(l)max(\,i) , 



uniformly in A G K. In fact, for \B\ < 1 this follows from the comparison with the integral using the 
Euler-MacLaurin formula 



v />oo / />oo \ 

£/(») = / /(^dx + O / \f"(x)\dx) 



and for |S| > 1 we can write the sum as J2k( l + B 2 (A/B - [A/B] - k) 2 )- 1 = 0(1 
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Using this and ( |4.6| ) we see that 

E {A - B^ h {m))- 2 <Cj2(Bm- B^ h (A/B))- 2 < Cmax (l, ±\ . 

m>\ m>l ^ ' 

Hence, uniformly in £ 2 , h 2 

V(A^ 1 (£ 2 / l2 CA,/ l2 (m 2 )-4^i6,h 1 (m 1 )))- 2 = 0(l)maxfl,-^) . (4.11) 

Proceeding as in the proof of Proposition || we introduce -F(A, £21 £>) as the number of (£1, hi, h 2 ) satisfying 
t-yh-y = £ 2 h 2 + p. We now have 

F(\,£ 2 ,p)<4 d(\£2h2+p\), 



e 1 h 2 + P ^o 



and 



00 

dof 



g(\,n)^ J2 E n\^p){ P )- 2 

p=-oo |^ 2 |< AT 
00 

E E d(\n\)d(\n+p\)(p)- 2 

p= — OO 1<M<«A 

(JVA \ 5 00 /N\+p \ 3 

J>(n) a E(P)" 2 E 
n=l / p=0 \ m=l / 

00 ] 

< C 3 (AiV(logA + logAf) 3 )=E (p)- 2 ((XN +p)log 3 {XN +p)y 

p=0 

< CAiV(logA + logiV) 3 . 
Using ( 4.11 ) we can estimate ( l.lOj ) by: 

C7max(l + |^|)| 5 (0|A- 3 J (C(A, A|£|) + l) < , 

completing the proof of the proposition. □ 

Remark. As was emphasized by the referee, Propositions || and || are essentially sharp. The "diagonal" 
solutions ji — j 2 , ki — k 2 , li — l 2 give a lower bound 1/A. Hence, no essential improvement of the Theorem 
in Sect.l is possible by this method. 

5. Construction of the family of surfaces 

To prove the main theorem we need to construct a family of surfaces, Q — {<? Q '' 3 }, for which in the 
expansion of the phase shifts (4.8), 

iP a 'P(x) = ax + P$(x), |$"t M _ e) |>C e >0, ( ,/3)€(ao-7,ao+7)x(-7.7), (5-1) 

so that we can appl y Propositions f| and ||. Recalling Proposition this means that we want to find a"^ 
satisfying (2.1)-(2.4), and such that 

r * {x) = I jf" ^1 _ ' _ (l _ x ^ dr = ax + /»(*) , (5.2) 

with $ convex or concave. 
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We will now skip the indices a and (3. If we write 



and 



then 



W(r) d ^^-l. 
a(r) 2 



(*) = / (5-3) 



W(r) 



7T \V1-X 2 J 



and we might study the simpler function (j> instead. From the assumptions on a, W is monotonically 
decreasing and r 2 W(r) is smooth and non-zero at r = 0. Hence there exists a smooth monotonically 
increasing function, y(r), such that 

Since we can write r as a function of y we define 

That way we can express 4>(x) as a linear transform of F: 

<f>(x) = xI(F)(x) , I(F)(x) = -lj F{xy)dy. (5.4) 

From this we immediately get a linear model correspoding to F(x) = t > 0: 



. . r 2 \J\ — x 2 , / x \ 1 



since 



((l-y- 2 )J-l)dy = -7r/2. 



Remark. The surfaces defined using the linear model do not have uniform pair correlations measures. In 
that case we can compute the leading contribution to E*(f) directly. To apply Propositions || and || we 
need to have the linear term in ip and that forces the singularity at for our surfaces: only one value of a 
corresponds to a smooth surface. 

We want to introduce the convex or concave term in ip by perturbing the case F = const. For that let 
us first establish some simple properties of the transform F i— ► 1(F). We denote by smooth functions 
with bounded derivatives and by Sp hg spaces of poly-homogeneous (classical) symbols. 

Lemma 2. For g e Cg°([0,oo)) 

I(g)(x) € Cg°([0,oo)) +xlogaCg°([0,c»)) . 

When g 6 5p h 2 g ([0, oo)) then 



I(9)\[i,oo)e S^([l,oo)) , I(g)(x) ~ (- f f(y)dy) - + £ |£ 

\ J / 1 n 



x — > oo . 

X K 

k=2 
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Proof. To prove the first part of the lemma we write 



I(9){x) = ((1 - y- 2 )\ - l) 9 (xy)dy + ((1 - y~ 2 )^ - l) g(xy)dy , C > 1 , 



o 

where the first term on the right hand side is clearly in C£°([0, oo)). In the second term, the integrand can 
be rewritten as 

Thus we are concerned with integrals of the form 

/ -g(xy)dy = x^ 1 Y~ k g(Y)dY + x^ 1 Y- k g(Y)dY . 

Jc y Ji Jcx 

where the first term is smooth and uniformly bounded in k. To study the second term we write 

9(Y) = go + giY + ■ --gi-iY 1 ' 1 + gi(Y)Y l 

which gives 

1 1 n-k+i 
' C~ k+1 9° + l 7;C- k+2 gix + ■■■ + r g l -ix l ~ 1 + x l F k ^ c {x) , 



k-1 k-2 k-l 

for k > I and 

r^-rC-^Jo + T^C- k+2 g lX + ■■■ + g k - X x k - x logx + x k - x G KW {x) , 
k — 1 k — 2 

for finitely many k < I. Since we check that 

G k ,i,c(x)eC 0o ([0,l/C)), k<l, 

\F k ,i, c (x)\ < G l (C- k+l+1 + \x\ k - l ~ x ) , k>l, 

we can sum up the contributions from different fc's (C > 1, |x| <C 1, and we use the uniform convergence of 
(1 — z) 5 ). Thus for every I we obtain 

I(g)(x) = h u (x) +xlogxh 2 ,i(x) + 0{x l ) , h hh h 2 ,i £ C°° , 

and consequently 1(g) <E Cg°([0,oo)) + xloga:Cg°([0, 00)). 

The second part of the lemma is even more clear. If for large Y, g(Y) ~ G k Y~ k , then then 

□ 

The lemma shows that we cannot expect smoothness of ip(x) at the end points x — 0,1 but that the 
function is very well behaved in the interior, as in any case is implicit in Proposition |^. 

Having discussed the general properties of the transform / we now state a straightforward 



Lemma 3. For = xI(f)(x/\/l — x 2 ) and < x(z) = z/y/l + z 2 < 1 we have 



&'(x(z)) = (1 + z 2 )^ 



((1 - ir 2 )l - l) ((2 + 3z 2 )yf(yz) + zy 2 (l + z 2 )f'{zy)) dy 



□ 

Guided by the two lemmas we can easily construct a family of surfaces for which tp has the needed 
properties. We want to find / e C£°(R) such that for j3 small enough and a close to ao, F a '^(y) = 
a + Pf(y) > 0, and so that a(r) obtained from inverting the process described above has the properties 
( 2.1 )-( 2.4). This is easily achieved by demanding that / is a symbol of order —2 on [0, 00). 
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We also want Q"(x) described in Lemma || to have a fixed sign. From the formula we see that that $ is 
concave if 

af'(y) + yf"(y)>0, y>0, a = 2,3. (5.5) 
In fact, the integrand has the same sign as 

(2 + 3z 2 )f'(Y) + (1 + z 2 )Yf"{Y)\ Y =y Z 



and that is positive for any z if (5.5) holds. 



We can summarize this discussion in 

a,0/ 



Proposition 6. For a in a neighbourhood of a fixed ao < and for (3 small enough, let a a ± (r) be obtained 
from f 6 "5^h g ([0j °°)) by the following procedure: 

F a >P(y)=-2a-[3f(y), 



dr 

Then a"' 13 has the properties (2A)-(2J:) and for the set of two parameter families of surfaces 



W^-(F a > f, U a > f, (r))- 1 , y Q ' /3 (0) = 0, a a ' (r) = y a < P {r){l + y a ' P (r) 2 )-^ 



G = {(X,gf ) : gf P = dr 2 + a a /(r)d9 2 , / g S^([0, oo)) satisfies (fcj), \a - a \ < S f , \[3\ < e f } 



the leading part of the phase shifts depends linearly on a and [3 and (5.1) holds. 

Combined with Propositions ^| and || this provides an infinite dimensional family of perturbations of the 
linear model each giving a two parameter family of surfaces for which the Theorem of Sect.l holds. Perhaps 
the simplest example is obtained by putting 

F°>< l {y) = -2a-p—}— 5 , < p < \ . 

1 + py z 6 

Proof of the Main Theorem. Proposition ^ guarantees that the leading parts of the expansions of the phase 
shifts of (X, g"' 13 ) satisfy the assumptions of Proposition |5[ Let p\ a ^ denote the pair correlation measure 
for (X, g a 'P) given by (L2). Recalling (4.1) and the discussion preceding it we see that all for / G S(R) 

PT'^U) = /(0) + /(0) + oa^oo(I) + E{(f; a, (3) . 

Proposition ^| now shows that the assumptions of Lemma [l] are satisfied with F(X) = log 3 A/A and that 
lemma gives the statement of the Main Theorem. □ 
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